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We review a new approach to calculating transverse momentum distributions of 
the Higgs and electroweak gauge bosons using the Soft-CoUinear Effective Theory. We 
derive a factorization theorem for transverse momentum distributions in terms of newly- 
defined impact-parameter beam functions (iBFs) and an inverse soft function (iSF). The 
iBFs correspond to completely unintegrated parton distribution functions and provide 
interesting probes of momentum distributions within nucleons. The numerical matching 
between the low and high transverse momentum regions is improved in this approach 
with respect to standard techniques. We present results for next-to-leading logarithmic 
resummation for the Higgs and Z-boson distributions and give a comparison with Teva- 
tron data. 



1. Introduction 

Low transverse momentum (pt) distributions of the Higgs and electroweak gauge 
bosons play an important role in numerous physics studies, including Higgs searches, 
precision measurements of the VF-boson mass, tests of perturbative QCD, and in 
probing transverse momentum dynamics in the nucleon. In the region px <C M, 
where M denotes the mass of the Higgs or electroweak boson in question, large 
double logarithms in p^/M spoil the convergence of the perturbative expansion 
in the strong coupling. A proper description of the transverse momentum distri- 
butions in this region requires resummation of these large logarithms. Extensive 
studies [Tl[lllll[3[Sl[71lili[Tni[ni[Tl[Tl[Tl[Tni[T5]oflowpT resummation have 
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been performed. Most are based on the well-known Collins-Soper-Sterman (CSS) 
formalism. 

More recently in Refs. [T71 [13 [H] , low pt resummation was studied using the 
Soft-CoUinear Effective Theory (SCET) [HI US- In the region M > pT > 
Aqcd^ the schematic form of the factorization theorem is 

^-H^g^f^f. (1) 

Here, _ff is a hard function that encodes the physics of the hard production vertex 
and whose renormalization group (RG) evolution sums the logarithms of p^/Af. Q 
is a gauge invariant Transverse Momentum Function (TMF) at the scale fix ^ Pt 
that describes the dynamics of initial state soft and coUinear emissions that recoil 
against the heavy boson. / denotes the standard PDF and is evaluated at /iy ^ Pt 
as determined by DGLAP evolution. In the region where pr ~ ^qcd, the TMF 
beomes non-perturbative and encodes the transverse momentum dynamics in the 
nucleons and of the soft radiation in the process. 

A smooth transition between the regions of non-perturbative and perturbative 
Pt is achieved by implementing a model ^ for the TMF Q in the non-perturbative 
region such that it has the correct RG evolution properties and smoothly reduces 
to the perturbative result for increasing pT . The field-theoretic definition of G pro- 
vides further insight into the structure of the factorization theorem. In particular, 
we identify new objects called Impact-parameter Beam Functions (iBFs) which cor- 
respond to fully unintegrated PDFs. The form of the factorization in terms of these 
iBFs has differences compared to the more traditional TMDPDF formalism (see 
other proceedings [531 [Ml 123 of this workshop for the most current overview). We 
give a brief discussion on this in the next section. In the following sections, we 
also give an overview of the factorization theorem, numerical results, and provide 
comparison with data. 



2. Overview of Factorization Theorem 

The first step in the derivation of the factorization theorem is to integrate out 
the hard production scale by matching QCD onto the SCET which describes the 
dynamics of initial state coUinear and soft emissions that recoil against the heavy 
boson. These modes in the SCET have the momentum scalings 

p„ ^ M{r]'^,l,T]), - M(l, 77^77), Ps ^ M{Tj,r],'q), V^^, (2) 

corresponding to the n-collinear, fi-collinear, and soft modes respectively. The heavy 
boson acquires finite transverse momentum by recoiling against these soft and 
coUinear modes which have transverse momenta of order px- Using the soft-collinear 
decoupling property of the SCET, the dynamics of the coUinear and soft emissions 
are separated into two zero-bin subtracted E32ZI iBFs {Bn, B^) and a soft function 
(S) respectively. The zero-bin subtraction is necessary to avoid double counting the 
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soft region. Such a subtraction also appears in the TMDPDF formahsm'^Sl. Using the 
equivalence of the zero-bin and soft subtractions one can rewrite the factorization 
theorem in terms of iBFs (Bn,Bn) defined without the soft zero-bin subtractions 
and an Inverse Soft Function (iSF). In the rest of this article, iBFs will refer to the 
Bn , Bn functions defined without the soft zero-bin subtraction. For perturbative 
values of the transverse momentum px 3> ^qcDi the iBFs are matched onto the 
standard PDFs after performing an operator product expansion in ^qcd/pt- The 
Wilson coefficients of the iBF to PDF matching combined with the iSF give the 
TMF function Q . We schematically summarize these steps in deriving the factoriza- 
tion theorem below: 

^^JpS\Mqcd\' (3) 
i (match QCD to SCETp^) 



dp^dY 



PS\C<E>{Oscet)\'' 

I (SCET soft-coUinear decoupling) 

- (g) B„ (g) (g) S" 

4, (zero-bin and soft subtraction equivalence) 

- (g B„ (g) Bfi (g) S*"^ 

i ( iBF to PDF matching) 

- ® [l„ g) Ift g) S-^] g)/, g) fj, 

Q 

where the coefhcients In,n arise from the iBF-to-PDF matching 

B„ = X„ ® /, Bri = Ifi (g /. (4) 



For pt ^ ^QCD, the iBFs, which correspond to fully unintegrated PDFs, describe 
the evolution of the initial state nucleons followed by their disintegration into a jet of 
collinear radiation that recoils against the heavy boson. Analogous beam functions 
are known to arise in other processes at hadron colliders [291 ED] ■ For pT ^ J^qcd , 
a perturbative matching of the iBF onto the PDF is no longer valid. In fact, in 
this region, the iBFs and the iSF are non-perturbative. Eq.Q should no longer 
be viewed as a perturbative matching but instead must be viewed as a definition 
of the new non-perturbative functions X„ and !„. Thus, the TMF function Q is 
non-perturbative for px ^ ^qcd- It encodes the transverse momentum dynamics 
within the nucleon and in the soft radiation that determines the px of the heavy 
boson. 

For the sake of brevity, we focus most discussion here specifically for the case 
of the Z-boson px distribution. We refer the reader for details of the derivation, 
notation, and the case of Higgs boson distribution to Refs. [TTIIIH]- The factorization 
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theorem for the Z-boson distribution takes the form 

. f ^4 f ^4 (5) 

dp^dr N^Jq Jo J^^ x[ J^^ 4 

X H'^{xiX2Q'^,fJ.Q;HT) g''''''{xi,X2,x[,X2,PT,Y,fJ,T)fr{x[,fJ,T)fs{x2^fJ.T), 

where the indices r, s run over the initial partons. The TMF function C/'J''* is given 

by 

/(Pb f 
j^;^jQ[bj_PT] J dt+dt^ 

X In;qr{^,tn,b±,fJ-T)In-qs{^,t^,b±,HT) (6) 
X^ X2 

tt 



The iSF ^^^^(wi, J)2, A^) is given by 

db+db 

where the position-space soft function Sqq(b^ ,b^ ,b±, fj,) is given by the vacuum 
matrix element of soft Wilson lines as 

Sgg{b+,b',b^,fi)^TT{0\f[SlSn]{b+,b-,b^)T[SlSnm\0). (8) 

The n-coUinear iBF is defined as a nucleon matrix element of collinear quark fields 
and Wilson lines in SCET as 

Bl{x,t,bj_,f,) = J '^e^^'-{p,\ir.qWAb-M^ 

X 5{n-V ^xn-p^)Wi£,^qmPi) (9) 

and the coefficients In-qr appearing in Eq.Q are given by the iBF to PDF matching 
equation 

dz ( x 

Bl{x,t,bj_,fl) = / — iln-qq'i-,t,bj_,^i) fg,{z,^i) 
Jx Z ^ Z 

X 1 

t,b^,^i) fg{z,iJL)\. (10) 

Analogous equations exist for the n-coUinear functions. 

As can be seen from Eq.([9]), the iBF corresponds to a fully unintegrated PDF. 
This is in contrast with the TMDPDF which has been extensively studied in the 
context of transverse momentum distributions. The TMDPDF is typically defined 
with reference to an external regulator to regulate spurious rapidity divergences 
that arise in perturbative calculations. The independence of the cross-section from 
this external regulator gives rise to the CoUins-Soper evolution equation which sums 
large logarithms. In contrast, the iBF is free of rapidity divergences. They are au- 
tomatically regulated by the additional residual light-cone momentum component 
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{k~^) which appears through the definition of the convolution variable t = Qk^. 
This is a reflection of the fact that for a finite pr gluon recoiling against the heavy 
boson, the mass-shell condition p^p'' = p^ ensures a non-zero residual light-cone 
momentum component so that no rapidity divergence occurs in the physical pro- 
cess. The use of the residual light-cone momentum component to regulate rapidity 
divergences differs from the TMDPDF approach where an external regulator is in- 
troduced instead. The iBF is more differential than the TMDPDF and corresponds 
to a fully unintegrated PDF. A more detailed understanding of the relationship 
between the TMDPDF and the iBF approaches requires further study. However, an 
explicit comparison of the iBF approach with the TMDPDF or CSS formalism is 
possible by expanding the resummed results and comparing the logarithms gener- 
ated at any fixed order in perturbation theory. Such an explicit check of the leading 
and next-to-leading logs in the iBF approach was performed in Ref . [18] . 

3. Numerical Results 

In this section, we present numerical results for the Higgs boson and Z-boson px 
distributions. We refer the reader to Refs. [T71 [TH] for analytic expressions for the 
iBFs and the iSF. In Fig.[l] we show the fixed order, Leading- Log (LL), and Next-to- 
Leading Log (NLL) resummed results for Higgs production at y/s — 7 TeV and zero 
rapidity for a Higgs mass of nih = 165 GeV. The plot is cut off at pT = 3 GeV so 
that the distribution is given entirely in terms of perturbatively calculable functions 
and the standard PDFs. Study of the distribution for non-perturbative values oipT 
requires a model for the TMF function. We describe our modeling procedure later 
for the case of Z production. From Fig.[l] we see that the effect of the resummation is 
dramatic and changes the shape of the distribution. In particular, the 1/px singular 
behavior of the fixed result is brought under control by resummation. 

Similarly, we present NLL resummation results for the Z-boson distribution for 
perturbative values of pr in Fig. [2] Also given is a comparison with Tevatron data 
collected by the CDF and DO collaborations [SI [32]. The plot in Fig. [2] is cutoff 
at pt ^ 1-75 GeV. As seen, good agreement is found with the Tevatron data. A 
description of the Z-boson px spectrum in the region px ^ ^qc d requires a model 
for the Z-boson TMF function. 

In Ref. [19], a phenomenological study of TMF models was carried out. The 
guiding principles for the construction of TMF models include reproducing the 
correct RG evolution properties and requiring that the TMF model reduces to the 
perturbative result as one increases px- These principles are encoded by writing the 
the TMF as a convolution [33ll34l|35l|36] of the perturbative expression of the TMF 
with a non-perturbative model function as 



G'^''''{xi,X2,x'-^,x'2,PT-,Y,llx) = / dp'xGps'rt.i.Xl,X2,x[,x'2.,PTjl + {p'x/PTy,Y,^J■T) 





X Gmod{p'x,0,^^)^ 



(11) 
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Fig. 1. Numerical predictions for the transverse momentum spectrum for Higgs boson produc- 
tion at the LHC for central rapidity. Shown are the fixed-order result and those obtained after 
implementing the resummation formula of Eq. ([sjl through LL and NLL. The bands arise from the 
scale variation. 



where Gmod denotes a model function that peaks near pip ~ Aqcd with parameters 
a, A that can be fit to data. The scale dependence of the TMF is contained entirely 
in the perturbative expression ^p^rt. of the TMF so that the required RG evolution 
properties are reproduced. For phenomenological analysis we parameterize Gmod as 

23/2-a 



G,„od(p^,a, A) = 



1 



A r(a-l/2) 



with the normalization condition 



poo 

/ dp'rp Gmodip'r, a, A) = 1 

JQ 



\ a— 1 


r p't'] 


1 exp 


. 2A2. 


= 1. 





(12) 



(13) 



A sensible choice for /i-r that can be applied in both the perturbative and non- 
perturbative pt regions is 



2 /■2 2 , 2 

= ? Pt+ Pt„ 



(14) 



where PTmin > 1 GeV is a low, but still perturbative, scale and can be viewed 
as another parameter of the model. ^ is a scale variation parameter we take to be 
0(1). The above choice of scale for /i^ has several useful properties. As px ^ 0, the 
scale Pt PTmin so that Gpart ^'i- (111 is ^till evaluated at a perturbative scale. 
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Fig. 2. Numerical predictions for the transverse momentum spectrum for Z boson production at 
Tevatron Run 1, compared with data from both CDF and DO. Shown is the resummation prediction 
at NLL accuracy. The bands arise from the scale variation, while the result for the central scale 
choice is shown by the solid line. The lower limit of the plot is PT= 1-75 GeV. 

Similarly, the running of the hard function H'^(xiX2Q'^ , ^J^q', fJ^x) will freeze at the 
perturbative scale pTmin as pr — 0. For larger values of pr ^ PTmm, l^T ^ £,Pt 
so that the appropriate choice of ~ Pt in the perturbative region is recovered. 
For Pt ^ ^qcd , the model TMF reduces to the expected perturbative expression 
up to power corrections as 

G'^''' {Xi, X2, x\, x'2,Pt,Y, fir) . = gpirt . ixi,X2,x[,x'2,PT,Y, (Xt) + 0{ ^^^° ), 



(15) 



so that a smooth transition between the non-perturbative and perturbative re- 
gions is achieved. In Fig. [3] we show the Z-boson pT spectrum including the non- 
perturbative Pt region. The values of the parameters a — 2.50, A — 1.43 GeV in 
Gmod used in Fig.jSj are obtained from a best fit to the CDF and DO data for pT < 10 
GeV. We see that the form of the model TMF is flexible enough to describe data 
over the entire range of pt- 



4. Conclusions 



We have presented a new factorization theorem for the low transverse momentum 
distributions of electroweak gauge and Higgs bosons using the SCET. The fac- 
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Fig. 3. The result for the py-spectrum of the Z-boson for the best fit parameter choices a = 
2.50, A = 1.43 GeV. We have varied /^t.Q within the range 1/2 < 5t,Q < 2 for fj,^ = £,tPx + 
PTmjn>/^Q = "SQ^I fol" PTrmn = 2 GeV. 



torization theorem is formulated in terms of Impact-parameter Beam Functions 
(IBFs) which correspond to fully unintegrated PDFs. We present results of NLL 
resummation for Higgs and Z-boson distributions. Good agreement is found with 
data collected by the CDF and DO collaborations at the Tevatron. More work is 
in progress to achieve NNLL results; recently in Ref. [37], NNLO results for the 
soft function of the Z-boson px distribution were obtained. We look forward to the 
further development of our formalism and applications to phenomenology. 
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